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Abstract 

Theories with large extra dimensions provide the possibility that a flavor neutrino, local- 
ized in a 3+1 brane, can mix with a singlet neutrino living in the bulk. This mixing leads 



to unconventional patterns of neutrino matter oscillations and we examine in details how 
these oscillations depend upon two parameters: the brane-bulk coupling £ and the effec- 
tive mass /i of the flavor neutrino inside matter. We find that high energy (E > 50 GeV) 
Up neutrinos, to be detected by neutrino telescopes, can give signals of extra dimensions. 
With a 1 km 3 neutrino telescope extra dimensions with radius down to Ifim can be tested 
directly, while for smaller radius an indirect evidence can be established. 



Introduction. The standard model (SM) of strong and electroweak interactions has 
been extremely successful. It provides a consistent theoretical framework within which 
we can analyze and understand all available experimental data. Still we know that it 
cannot be regarded as the final theory. Any attempt to include quantum gravity leads to 
a unified theory where two disparate scales coexist, the electroweak scale (Myy ~ 1 TeV) 
and the Planck scale (Mpi ~ 10 19 GeV). Quantum radiative corrections then, especially 
in the scalar sector (Higgs field) of the theory, tend to mix the scales and without an 
incredible amount of fine tuning, will always equalize them (hierarchy problem). 

A novel approach has been suggested to alleviate the hierarchy problem IJ. Our 
four-dimensional world is embedded in a higher dimensional space with D dimensions 
(D = 4 + n). While the SM fields are constrained to live on the 4-dimensional brane, 
gravity can freely propagate in the higher- dimensional space (bulk). The fundamental 
scale Mf of gravity in D dimensions is related to the observed 4-dimensional Planck scale 
Mpi by 

M 2 Pl = Mf n V n (1) 

where V n is the volume of the extra space. For a torus configuration 

V n = $w) n R 1 R 2 ...R n (2) 

with Ri (i = 1, 2, ...n) the radii of extra dimensions. Then for a sufficiently large volume 
V n the fundamental scale of gravity Mf can become as low as M w . In this radical way 
the hierarchy problem ceases to exist as such. 

Besides the graviton, a sterile neutrino can freely propagate in the extra dimensions. 
Since the extra dimensions are compact, a bulk neutrino will appear at the 4-dimensional 
brane as a Kaluza-Klein (KK) tower, i.e. an infinite number of 4-dimensional spinors. The 
Yukawa coupling of the standard lepton doublet, the Higgs scalar and the right-handed 
bulk neutrino will provide a mixing between the left handed neutrino of the standard 
model and the KK modes. In the next part we present this mixing and we obtain the 
eigenvectors and the eigenvalues of the mass matrix. In the third part we focus our 
attention into neutrino oscillations, inside matter. Compared to the usual oscillations, 
novel features appear since now we have a coupled system of infinite degrees of freedom. 
The pattern of oscillations depends upon two parameters : the coupling £ between the 
left-handed neutrino and the KK states and the effective mass fi of the flavor neutrino 
inside matter. We study in detail this dependence. In the final part we suggest that 
high energy cosmic neutrinos can probe more efficiently the bulk of the space. These 
neutrinos, after traversing the Earth, will be detected at future neutrino telescopes. We 
present detailed spectra of the comic as function of the energy and the incident angle. 
At the end we present our conclusions . 



The Brane-Bulk Model. The available gravity experiments imply that any devia- 
tion from the Newton law should occur at distances r less than 1 mm. Then, to satisfy 
equations (l)-(2), we need more than two extra dimensions . To simplify our investigation 
we consider a single large extra dimension, assuming all other dimensions smaller and not 
affecting our analysis. 

There is already an extensive literature on neutrinos living in extra dimensions [0]-|[3 



We start by considering the action for a massless bulk fermion \1/ living in 5 dimensions 

S 9 = J d A xdy^T»d^ M — 0,1, ...5 (3) 

\I> is decomposed as 



N L 

each component admitting a Fourier expansion 



*=(^ R ) (4) 



N{x, y) = -^=N (x) + jrJ=( Nn (x) cosffi + JV„(x) sinffl) (5) 

n=l x ' 

Thus at the brane the bulk fermion appears as a KK tower. Since the momentum in the 
exrtra compact dimension is quantized, each KK mode appears as having a mass n/R 
(n = 1,2, ...). Consider now the coupling of the standard left-handed lepton doublet to 
the bulk neutrino 



U LHN R 5(y) (6) 



where H is the Higgs scalar doublet and h a dimensionless Yukawa coupling. After the 
Higgs field develops a vacuum expectation value (v ) , we get the mass terms 



mb>L 



(N m + V2Y,NRn) (7) 



n=l 



with m ~ hvM f /Mpi. Notice that for M f ~ ITeV 7 , h ~ 1, we obtain m ~ 10~ 4 eV 0-0. 
Also, the left handed zero mode A^^ decouples from the spectrum and remains massless. 
Denoting 



the mass term is ^ r with the mass matrix M given by 



M 



( m \/2m y/2m 

1/R 

2/R 

V ' 



(9) 



The evolution of the neutrino states is determined by the Hamiltonian 



H 



2E U 



-MM 1 



(10) 



In terms of rescaled variables, H = X/(2E U R 2 ), where the dimensionless matrix X be- 
comes (we keep k rows and columns) 



X 



( (2k + 1)£ 2 v^f 2-v/2£ 3\/2£ 

V2£, 1 

2y/2£ 4 

3y/2£ 9 



V 



(k-i)y/2£ 



(*-l)V2£\ 









(Jfe-1) 



/ 



with £ = mi?. For a neutrino propagating inside matter the upper-left element is replaced 
(2k + 1)£ 2 — ► (2/c + 1)£ 2 + fi, where fi is proportional to the effective mass of the neutrino 
inside matter. For a neutrino [i = \f2E v R 2 GpN n with N n the neutron density. The 
eigenvalues A 2 of the matrix X satisfy the equation 



At - X +£ (A7r) cot(A7r) 



n^-A 2 )=o 



(12) 



71=1 



The corresponding eigenvector to A 2 is £>„ = (e n0 , e n i, e n2 , •••) with 



(13) 



e n0 I - + 



2A 2 



+ • 



(A n - ^ 2 



A n ' 



2? 



(14) 
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Neutrino Mixing. Before presenting the numerical solution to equation (|T|) it is 
useful to mention some qualitative features 

i) small £ 

For £ = there is no coupling between the SM neutrino and the KK states. The left- 
handed neutrino remains massless, Ao = 0, while for a neutrino propagation inside matter 
■^o — y/V- The other eigenvalues are \ n = n (n = 1,2,3,...). As £ starts to rise 
A n ~ n + £ 2 / (nn) (for large n) 

ii) large £ 

For large £, A n tend to n + 1/2. An expansion for large £ and large n provides 



X n^(n + -)(!-_) (15) 

iii) resonance structure 
The amount of admixture of the n-th eigenstate B n to the neutrino state v is determined 
by |e„o| 2 and it is becoming maximal when the resonance condition 

/in = A 2 (16) 

holds (see equation |l4j) . Since A„ ~ n, the resonance condition is simplified to p ~ 
n 2 (n = 1, 2, 3, ...). For small £ values the resonance is narrow, while for large £ values 
the resonance is broad and more than one state might be involved. Recall that for a 
neutral medium 



_ G F R 2 pE v 
^~ V2M N 

with p the density of the medium and the nucleon mass. The requirement p ~ 1 
(strong oscillation phenomena) leads to a relationship between the model parameters (R), 
the medium properties (p) and the neutrino energy (E u ) 

, P ^ , E» ^ , R ,2 ^ 1 m 



v 10 v 100 GeV v l am' 

cm" r~ 

iv) Level repulsion 

When there is no coupling among the SM neutrino and the KK states (£ = 0), the eigen- 
values cross each other (see dashed lines in figure [I] ). When the off diagonal elements 
become non-zero (£ 7^ 0), the well-known phenomenon of level repulsion appears : as the 
lower eigenvalue increases by increasing p, by approaching the resonance region (// ~ 1) 
levels off and reaches asymptotically the value 1, thus avoiding crossing the next eigen- 
value. Similarly the higher eigenvalue starts at A ~ 1, increases as p increases and when 
encounters the second resonance (p ~ 4) levels off at the value 2. (see the figure [I] ) 

v) graphical solution 

The eigenvalues A 2 can be found graphically as the intersection points of the two curves 
(see the figure || ) 



y 1 = A 2 - n 



(19) 



y 2 = £ 2 (\tt) COt(ATT) (20) 

Most of the qualitative figures described above can been seen by inspecting the figure 



Oscillating Neutrinos and Neutrino Telescopes. From the previous analysis it 
is clear that for neutrinos traveling inside the Earth, strong oscillation phenomena occur 
whenever the condition ( jl8f) is satisfied. Thus for R ~ 1/im, we need neutrinos of energy 
E ~ 100 GeV, while a smaller radius is probed by even higher energy neutrinos. It 
is clear that neutrino telescopes |14| which are able to detect neutrinos from ~ 20 



GeV up to ~ 50 TeV are best suited for this study. At this early stage of investigation 
we adopt a simple model where the flavor state mixes only with a bulk neutrino. We 
expect this simple model will indicate the salient features of the more complicated physical 
phenomenon. 

Once a flavor state is produced, its time evolution will be governed by the equation 
W)>= <oexp(;^|) \B n > (21) 

n 

where the summation over n is carried up to n max = [E V R] The amplitude for the transi- 
tion Up — > is 

A=Y, I e -I 2 exp(^|) (22) 



Figure § shows the probability Piy^ — > v^) for a neutrino traversing a medium of constant 
density, as a function of distance. The value we used for the density of the medium 
corresponds to the density of the Earth's mantle. We used also R = lfim. A clear dip 
develops at the appropriate length scale. The probability Piv^ — > v^) as a function of 
energy is presented in figure |] . Notice that an impressive resonance structure appears at 
the energy which corresponds to \i ~ 1 and all disappear. At the resonance, reading 
the value of the energy and using equation ( p~8|) we can extract the value of the radius R. 

It is useful to evaluate the transition probability for a neutrino traversing the Earth. 
For the Earth's density we adopt a model where the core is a sphere with constant density 
11 gr/cm 3 while the mantle, surrounding the core has constant density 4.4 gr/cm 3 ||T5| . 



We consider a neutrino which for time ti propagates inside the mantle, for time £2 crosses 
the core and it takes time £3 within the mantle again . The amplitude Aiv^ z/ M ) is 
given by 



n,l,r 

where tilded (untilded) parameters correspond to core (mantle) and D is the rotation 
matrix from one set of eigenvectors to the other 



D nl = e ^*ik (24) 

k 

The probability Piy^ — > v^) for a neutrino going through the center of the Earth, as 
a function of the energy is shown in figure [5|. The resonance dip is now displaced. 

How efficiently can a neutrino telescope unravel the underlying dynamics? Without 
entering into technical details and specifications of a neutrino telescope, we can recall 
that by determining the direction of the produced muon, we also fix the direction of the 
parent muon neutrino. Thus we can have a zenith coverage of incident neutrino from 
9 = (neutrinos going through the center of Earth) to 9 ~ tt/2 (neutrinos just scratching 
the Earth, before being detected). In the absence of any oscillations the neutrino flux is 
^-independent. If oscillations are present, by varying 9 we vary the distance x traveled 
by i/p (x varies continuously from 10 km for 9 = tt/2 to 13000 km for 9 = 0). In our case 
the 9 ~ tt/2 neutrino flux remains undistorted by matter effects and represents therefore 
the initial incident neutrino flux. The deduced ratio I (9) /I (9 = tt/2) (1(9) is the infered 
intensity of neutrino flux at angle 9) reflects then the probability for a muon neutrino to 



oscillate. The energy of the incident neutrino is determined for the contained events \W 
Thus we have a set of events, where both the distance x and the energy of the neutrino 
E v are determined, allowing to extract the theoretical parameters (£ and R). Neutrino 
induced events with energy up to 300 — 400 GeV are expected to be contained within a 
1 km 3 volume and therefore we deduce that 1 km 3 neutrino telescope can provide direct 
evidence for an extra dimension with R > lfim. For smaller R, we expect the resonance 
structure to shift to higher energies and thus the resonance dip will not be directly seen. 
Still we may have an indirect evidence for strong oscillations by examining the ratio 
r = E s /E p , with E s (E p ) the number of stopping (passing) muons |17]]. If at 9 = the 



number of highly energetic (E v > 400 GeV) is depleted, then we expect the ratio r to 
increase for 9 = compared to the same ratio for 9 ~ tt/2. 

Another important scale in the neutrino oscillations is the resonance oscillation length 
L R . A significant depletion of occurs whenever the distance traveled x is equal to the 
oscillation length L, defined by 

L = 2n^- 2 (25) 



At the resonance , \x ~ 1, L becomes 

2^2 it M N 



■iR 



G F P (A? - Ao) 



(26) 



Whenever both equations fll!j|) and (^) are satisfied, we observe the maximum deple- 
tion of {P{Vfx — > v^) — 0). The eigenvalues difference Af — Aq depends upon £. For 
small £ values the resonance is narrow and the difference A^ — Aq is small. The corre- 
sponding Lr is big, much bigger than the Earth's diameter. For example for £ = 0.05 and 
with p = 10gr/cm 3 we find Lr = 23571 km. In that case strong oscillations cannot be 
observed. For £ = 2.2 (the value we opted), we find Lr = 6000 km, smaller than Earth's 
size, and strong oscillations are visible. 

Conclusions. We analyzed the situation where the left-handed neutrino z/ M of the 
SM mixes with a right handed neutrino experiencing a large extra space dimension. This 
mixing induces neutrino oscillations, even if the neutrino itself is massless. We encoun- 
tered novel features compared to the standard oscillations, since the mixing involves an 
infinity of KK states. A rich resonance structure appears accompanied by many spikes. 
We presented in detail the first resonance (/i ~ 1). A similar study can be realized for the 
second resonance (/i ~ 4). We emphasized that 1 km 3 neutrino telescope may be used 
as an instrument to explore the bulk space and establish unambiguously the existence 
of extra dimensions with R > Ifim. We didn't address the issues related to low energy 
neutrino data (solar neutrinos, atmospheric neutrinos), since for the R values we used, 
there is no significant mixing of low energy neutrinos with the bulk states. After all, the 
recent neutrino data (SNO |18[) disfavor the presence of a sterile neutrino. Altogether 
our analysis shows that neutrino mixing in models with large extra dimensions is an in- 
teresting phenomenon and deserves further study. 
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Figure 1: Eigenvalues A as a function of /i. Dashed lines for £ = (no coupling among 
the SM neutrino and the KK states) and the solid lines for £ = 0.2 
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Figure 3: The probability P{v^ — > z/ M ) for a neutrino of energy lOOGeV traversing a 
medium with density p = 4.4 gr/cm 3 (close to the density of the Earth's mantle) as a 
function of distance. 
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Figure 4: The probability Piv^ — > i^) for a neutrino traversing 10649 Km through the 
mantle of Earth ( with density p = 4.4 gr/cm 3 ), as a function of energy. 
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Figure 5: The probability Piy^ — > v^) for a neutrino going through the center of the 
Earth , as a function of energy. 
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